Measurement of double-spin asymmetries associated with deeply virtual Compton 
scattering on a transversely polarized hydrogen target 
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Abstract 

Double-spin asymmetries in exclusive electroproduction of real photons from a transversely polarized hydrogen target are measured 
with respect to the product of target polarization with beam helicity and beam charge, and with respect to the product of target po- 
larization with beam helicity alone. The asymmetries arise from the deeply virtual Compton scattering process and its interference 
with the Bethe-Heitler process. They are related to the real part of the same combination of Compton form factors as that determin- 
ing the previously published transverse target single-spin asymmetries through the imaginary part. The results for the double-spin 
asymmetries are found to be compatible with zero within the uncertainties of the measurement, and are not incompatible with the 
predictions of the only available GPD-based calculation. 

Key words: HERMES experiments, GPDs, DVCS, transversely polarized hydrogen target 
PACS: 13.60.-r,24.85.+p, 13.60.Fz, 14.20.Dh 



Preprint submitted to PLB 



September 16, 2011 



1. Introduction 

Generalized Parton Distributions (GPDs) provide a frame- 
work for describing the multidimensional structure of the nu- 
cleon HI IH 0] . They encompass information on the correlated 
transverse spatial and longitudinal momentum distributions of 
partons in the nucleon |4j, |5l |a |7|, [fy, |9[] . Furthermore, access to 
the parton total angular momentum contribution to the nucleon 
spin may be provided by GPDs through the Ji relation H. 

Hard exclusive leptoproduction of a meson or photon, leav- 
ing only an intact nucleon in the final state, can be described 
in terms of GPDs. GPDs depend on four kinematic variables: 
t, x, and Q 2 . The Mandelstam variable t = (p - p') 2 is the 
squared four-momentum transfer to the target nucleon, with p 
(p') its initial (final) four-momentum. In the 'infinite'-target- 
momentum frame, x and £ are related to the longitudinal mo- 
mentum of the struck parton, as a fraction of the target mo- 
mentum. The variable x is the average of the initial and final 
momentum fractions carried by the parton, and the variable 
known as the skewness, is half of their difference. The evolu- 
tion of GPDs with Q 2 = -q 2 , where q — k — k' is the difference 
between the four-momenta of the incident and scattered lep- 
tons, can be calculated in the context of perturbative quantum 
chromodynamics as in the case of parton distribution functions. 
This evolution has been evaluated to leading order JH B Bl [Toll 
and next-to-leading order dill El in the strong coupling 
constant a s . The skewness £ can be related to the Bjorken scal- 
ing variable .yb = Q 2 /(2p ■ q) through £ Xb/(2 - x B ) in the 
generalized Bjorken limit of large Q 2 , and fixed x% and t. There 
is currently no consensus about how to define £ in terms of ex- 
perimental observables; hence the experimental results are typ- 
ically reported as projections in xg. The entire x dependences 
of GPDs are generally not yet experimentally accessible, an ex- 
ception being the trajectory x = £, Il4ll5ll . 

The description of a spin- 1/2 hadron such as the nucleon in- 
cludes four leading-twist quark-chirality-conserving GPDs H, 
E, H, and E HSllEil. The GPDs H and E are quark-helicity 
averaged, whereas H and E are related to quark-helicity differ- 
ences. The GPDs H and H conserve nucleon helicity, while E 
and E are associated with a helicity flip of the nucleon. GPDs 
can be constrained through measurements of cross sections and 
asymmetries in exclusive processes such as exclusive photon 
or meson production. In the case of photon production, the 
asymmetries arise from the Deeply Virtual Compton Scattering 
(DVCS) process, i.e., the hard exclusive leptoproduction of a 
real photon, where the photon is radiated by the struck quark, 
and its interference with the Bethe-Heitler (BH) process, where 
the photon is radiated by the initial- or final-state lepton. The 
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DVCS process is currently the simplest experimentally accessi- 
ble hard exclusive process that can be used to constrain GPDs, 
but only on the trajectory x — £, . 

A variety of results from DVCS measurements at DESY 
and Jefferson Laboratory was published in the last few years. 
This includes results on beam-helicity and beam-c harg e asym- 



metries from CLAS ll7l LLSL ll9D and HERMES |gd 



21, 22] 



as well as cross-section measurements in Hall-A 12311 . all of 
which can serve to constrain mainly GPD H. HERMES addi- 
tionally obtained results on transverse-target asymmetries [24], 
which can constrain GPD E. Knowledge on both H and E 
opens access towards the determination of the total u and d- 
quark angular momentum through the Ji sum rule |3J]. This 
paper reports the first measurement of azimuthal asymmetries 
with respect to target polarization combined with beam helicity 
and beam charge, and with respect to target polarization com- 
bined with beam helicity alone, for exclusive electroproduction 
of real photons from a transversely polarized hydrogen target. 
One of these new asymmetries also has the potential in principle 
to constrain GPD E. 

2. Deeply virtual Compton scattering 

2.1. Scattering amplitudes 

The five-fold differential cross section for the leptoproduc- 
tion of real photons from a transversely polarized hydrogen tar- 
get reads IU6I1 



dV 



x B e" 



\T\ 2 



dx B dQ 2 d\t\ d0 d<f> s 32 (2nf yTT? 



(1) 



Here, e is the elementary charge, s = 2x B M p / y/Q 2 , where M p 
is the mass of the proton, and T is the reaction amplitude. Two 
azimuthal angles <p and <ps appear in the cross section in the 
case of transverse polarization of the target, and are defined in 
Fig.Q] 

The initial and final states of the DVCS process are indistin- 
guishable from those of the BH process. Hence the cross sec- 
tion contains the coherent superposition of the BH and DVCS 
amplitudes: 



1*7"! - \7~bh + TdvcsI 

= |7~bh| 2 + I^dvcsI 2 + 7~dvcs 7~bh + ^"dvcs ^"bh . 



(2) 



where T denotes the BH-DVCS interference term. The BH 
amplitude is calculable to leading order in QED using the form 
factors measured in elastic scattering. The interference term 
I and the squared DVCS amplitude |7~dvcsI 2 m Eq. |2] provide 
experimental access to the (complex) DVCS amplitude through 
measurements of various cross-section asymmetries II 1 611 . 

Each of the three terms of Eq. [2] can be decomposed as a 
Fourier series: 
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Figure 1 : Momenta and azimuthal angles for exclusive real-photon electropro- 
duction in the target rest frame. The quantity <f> denotes the angle between the 
lepton scattering plane containing the three-momenta k and k! of the incoming 
and outgoing lepton and the photon production plane correspondingly defined 
by the vector <f = k -k' and the momentum if' of the real photon. The symbol 
ips denotes the angle between the lepton scattering plane and S± , the compo- 
nent of the target polarization vector that is orthogonal to if. These definitions 
are consistent with the Trento conventions (25h and differ from those used in 
Ref. Q(J: = 7r- <p\^ and <p - (ps = n + 
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The symbols tf BH = ? (1 + e 2 ) ], ^dvcs = 1/(G 2 ) and 

/Ti = 1 /(xb y denote kinematic factors, where y = (p-q)f(p-k) 
and ei stands for the (signed) lepton charge in units of the el- 
ementary charge. Also, A — ±1 and S ± are respectively the 
helicity of the lepton beam and the magnitude of the vector 
the component of the target polarization vector that is orthog- 
onal to q. The BH coefficients c™ np , c™ p and s™ p in Eq. [3] 
depend on electromagnetic elastic form factors of the target, 
while the DVCS (interference) coefficients c°^ p s (cj, unp ), s°^ p s 

(4,un P ). <Sfp + V) (<tp+h) and 0»,tp + (-)) involve var- 

ious GPDs. The squared BH and interference terms in Eqs. [3] 

and |5] have an additional <p dependence in the denominator due 
to the lepton propagators V\ ((f)) and Vij^i) | 3 2^1 ■ The Fourier 
coefficients appearing in the interference term can be expressed 
as linear combinations of Compton Form Factors (CFFs) flill . 
while the coefficients from the squared DVCS term are bilin- 
ear in the CFFs. Such CFFs are convolutions of corresponding 
GPDs with hard scattering coefficient functions. 

The coefficients of particular interest in this paper are c™ p , 
s?5 D , c D Y,!r s , 5°v^ s c 1 and s 1 ^ u , which relate to double- 

l,lr' /i,TP+ 1,TP— ' n,ir+ n,lr— ' 

spin asymmetries involving transverse target polarization. The 
subscript 'TP' is used for BH terms, while the subscript 'TP+' 
('TP-') is used for DVCS and interference terms containing 
cos(0 - <ps) (sin(0 - (ps))- (The dependences of beam-charge 
and charge-difference or charge-averaged single-spin asymme- 
try amplitudes on remaining Fourier coefficients in Eqs. [3][5] 
were discussed in previously published HERMES papers B2.CX 
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Here, A" is a kinematic factor and the C-functions C l TP+ , C\ p _ 
and AC^ p+ can be expressed as linear combination of four CFFs 

( ( H, S, ( H, and <5) and the Dirac and Pauli electromagnetic form 
factors F\ and F 2 : 
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Note that even if the cross sections were measured for all eight 
possible combinations of beam charge and helicity and target 
polarization, at fixed x B and Q 2 it would be impossible to sep- 
arate the coefficients c™ p (s B j P ) and c°^,pf (s°jp S ). Neverthe- 
less, the BH coefficients can be calculated from the measured 
elastic form factors. 



2.2. Azimuthal asymmetries 

The asymmetries in the cross section for scattering of a lon- 
gitudinally polarized electron/positron beam off a transversely 
polarized hydrogen target, which embody the essential features 
of the Fourier coefficients appearing in Eqs.[3]|5] can be defined 
through 



do- = dcr uu (0)|l + e c m c ((f>) + ^°u CS (0) + S x & 
+ e € ml v (cf>) + e c S ± 3l\j T (<f>, cf> s ) 
+ AS x^ +DVCS (0, <f> s ) + e e AS x Ml T (<f>, </> s ) 
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where do-uu is the cross section for an unpolarized target av- 
eraged over both beam charges and both beam helicities. Us- 
ing the cross sections defined for purely polarized target states 
QS ±\ - 1) these asymmetries are expressed as 

&c(<f>) = — l - — \(da +1 + da +l + da +1 + do* 1 ) 
8do- U i] L 

-{da^ + do-- [ + do r ^ +do r - [ )\, (13) 



oda uu 

+ {da^ + do-' 1 - da' 1 - da' 1 )] , (14) 
odo- uu 

- (da' 1 + da~ l - do* 1 - do 7 ' 1 )] , (15) 



= — !— [(<^ + t _ j#H + _ d^+i) 
8do- vv L 

+ {do-- 1 - do-- [ + d^ 1 - da' 1 )] , 
?—\(d2 +1 - da H + d<F +1 - da H ) 



(16) 



- (da' 1 - do-' 1 + da' 1 - d^ 1 )] , 



(17) 



(10) M 



BH+DVCS 



( 0, s) = _^_[ ( ^+t _ j^n _ d ^ + d pn } 

8do- V ij L 

+ (do-- 1 - do-' 1 - do* 1 + da- 1 )] , (18) 



Mi T ((/>,4>s) = 7r-j~- — \(da +1 - do- +i - do* 1 + da +i ) 
odo- uu 

- (da' 1 - da' 1 - d^ 1 + do 7 ' 1 )] , (19) 

where the symbol + (-) denotes positive (negative) beam charge, 
— » (<— ) positive (negative) beam helicity, and f (J,) the target 
transverse-polarization direction. The arguments <f> and 0s are 
suppressed on the right-hand sides for brevity. 

3. The HERMES experiment and event selection 

The data reported here were collected with the HERMES 
spectrometer B27I1 using a longitudinally polarized positron or 
electron beam of energy 27.6 GeV scattered off a transversely 
polarized hydrogen gas target internal to the HERA lepton stor- 
age ring at DESY. The lepton beam was transversely self po- 
larized by the emission of synchrotron radiation 12811 . Lon- 
gitudinal polarization of the beam at the target was achieved 
by a pair of spin rotators in front of and behind the experi- 
ment 112911 . The sign of the beam polarization was reversed ap- 
proximately every two months. Two Compton backscattering 
polarimeters 1 3(3, 31 ] measured independently the longitudinal 
and transverse beam polarizations. 

The target cell was filled with nuclear-polarized atoms from 
an atomic beam source based on Stern-Gerlach separation with 
radio-frequency hyperfine transitions 1I32I1 . The nuclear polar- 
ization of the atoms was flipped on a time period of 1-3 min- 
utes. The polarization and the atomic fraction of the target gas 
were continuously monitored (33l 13 Si ■ The average values 
of the longitudinal beam polarization P{ and transverse target 
polarization S j for the various running periods are given in Ta- 
ble Q] Beam polarization and luminosity are given for the two 
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Lepton 


Longitudinal Beam 


Transverse Target 


Luminosity 


type 


Polarization (P ( ) 


Polarization (S t) 


[pb- 1 ] 












e~ 


- 0.286 


+ 0.338 


+ 0.721 


29.1 20.1 


e~ 


-0.286 


+ 0.338 


-0.721 


28.9 20.6 


e + 


-0.401 


+ 0.323 


+ 0.721 


11.8 17.5 


e + 


-0.401 


+ 0.323 


-0.721 


11.7 17.6 


Total 








81.5 75.8 



Table 1 : The type of the beam particle, the luminosity-averaged beam and target 
polarizations and the integrated luminosity of the data sets used for the extrac- 
tion of the various asymmetry amplitudes on a transversely polarized hydrogen 
target. The data were taken with an e + beam during the years 2003 (6.9 pb" 1 ) 
and 2004 (51.7 pb" 1 ) and an e" beam during 2005 (98.7 pb" 1 ). The uncertain- 
ties for the beam and target polarizations are 2.2% and 8.3%, respectively. 



beam-helicity states separately. The statistical uncertainties of 
the results reported here are generally larger than those reported 
in Ref. 12411 because here they scale as the inverse of the beam 
polarization. The target-polarization component S j is orthogo- 
nal to the direction of the incident lepton beam, while S j_ is or- 
thogonal to the direction of the exchanged virtual photon. This 
distinction is neglected in this analysis. 

The scattered leptons and produced particles were detected 
by the spectrometer in the polar-angle range 0.04 rad < 9 < 
0.22 rad. The average lepton-identification efficiency was at 
least 98% with hadron contamination of less than 1%. 

In this analysis, it was required that events contain exactly 
one charged-particle track identified as a lepton with the same 
charge as the beam lepton, and one photon producing an en- 



ergy deposition E y > 5 GeV in the calorimeter and > 1 MeV 
in the preshower detector. The following kinematic require- 
ments were imposed on the events, as calculated from the four- 
momenta of the incoming and outgoing lepton: 1 GeV 2 < Q 2 < 
10 GeV 2 , W 2 > 9 GeV 2 , v < 22 GeV and 0.03 < x B < 0.35, 
where v = (p ■ q)/M p and W 2 = M 2 + 2M p v - Q 2 . The angle 
between the laboratory three-momenta q and q' was limited to 
be less than 45 mrad, and —t < 0.07. 

An 'exclusive' event sample was selected by requiring the 

2 , with p = (M p , 0, 0, 0). As the data 
sample analyzed here is contained in that used in Ref. l24ll . 
missing only the 7.5% of that data set recorded in 2002, the M 2 , 



squared missing mass M 2 = (q + p - q') to be close to the 
squared proton mass Mp, with p = (M p , 0, 0, 0) 



distribution is very similar to Fig. 3 of Ref. II24I1 . The 'exclu- 
sive region' for e + data is chosen to be -(1.5) 2 GeV 2 < M^ < 
(1.7) 2 GeV 2 S. This region was shifted by 0.18 GeV 2 for the 
exclusive events from e~ data. This shift corresponds to the ob- 
served difference between the M 2 , distributions of the e~ and e + 



data samples 124 1. 



4. Extraction formalism 

The simultaneous extraction of Fourier amplitudes of beam- 
charge and target-spin asymmetries combining data collected 
during various running periods for both beam charges and helic- 
ities on a transversely polarized hydrogen target is described in 
Ref. [ 24] . It is based on the maximum likelihood technique 113711 , 



which provides a bin-free fit in the azimuthal angles <p and <f>s ■ 
In this paper, almost the same data set is analyzed, omitting the 
running periods when the beam polarization was small. This 
analysis differs in that the double-spin asymmetry amplitudes 
related to the ,?Ilt terms of the cross section given in Eq. Q~2] 
are also extracted. Furthermore, eight event weights were em- 
ployed in the fit to account for luminosity imbalances with re- 
spect to beam charge and beam and target polarizations, a tech- 
nique introduced in Ref. i24ll . 

Based on Eq. Q~2] the distribution in the expectation value 
of the yield for scattering of a longitudinally polarized elec- 
tron/positron beam from a transversely polarized hydrogen tar- 
get is given by 

{N){e c , P f ,S T ,<p,<f>s) = -C (e e , Pe, S T ) n(</>, <f>s) do- ut j(0) 



+ S r ^ cs (<^) 



x{l + e f ^c(0) + P^ 

,4>s) + ecPcS T 3ll T ((f>,<Ps)}, (20) 



+ eePe&lvW) + e ( S T W l m ((f>, 4>s) 



+ P e S T ^ +DVCS 



where £, denotes the integrated luminosity and 77 the detection 
efficiency. The asymmetries 3l c , ■^ J cs , ^°j CS , ^U' ^ut> 
^lt +DVCS ' anc ^ ^lt ^ related to the Fourier coefficients in 
Eqs.|3]-|5]and are expanded in terms of the same harmonics in <p 
and <p - 4>s m order to extract azimuthal asymmetry amplitudes 
in a maximum likelihood fit: 
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where the approximation is due to the truncation of the Fourier 
series. 

The amplitudes of beam-charge, beam-helicity and target 
single-spin asymmetries extracted in this analysis with 27 pa- 
rameters in the fit were compared with analogous results ob- 
tained with fewer parameters using the same Monte Carlo data 
sample. It was found that they agree with high accuracy. 

5. Background corrections and systematic uncertainties 

The asymmetry amplitudes are corrected for background 
contributions from the decays to two photons of semi-inclusive 
neutral mesons (mainly pions) and of exclusive neutral pions, 
using the method described in detail in Ref. 12411 . After ap- 
plying this correction, the resulting asymmetry amplitudes are 
expected to originate from single-photon production leaving the 
target proton intact as well as the associated production involv- 
ing excitation of the target proton (see the bottom row in the 
figures in the result section). Due to the limited resolution in 
missing mass and without detection of the recoil proton, the 
contribution of the latter process that falls within the exclusive 
window remains part of the measured signal. 

As the target polarization is involved in all the asymme- 
tries reported here and it was flipped on a time period of 1-3 
minutes, the effects of any time dependence of detector effi- 
ciencies or acceptance can be safely neglected. The dominant 
contributions to the total systematic uncertainty are the effects 
of the limited spectrometer acceptance and from the finite bin 
widths used for the final presentation of the results. The lat- 
ter is determined as the difference of the asymmetry amplitudes 
evaluated from yields integrated over one bin in all kinematic 
variables, compared to the amplitudes calculated at the average 
values of the kinematic variables. The combined contribution 
to the systematic uncertainty from limited spectrometer accep- 
tance, detector smearing, finite bin width, and imperfections 
in the alignment of the spectrometer elements with respect to 
the beam is determined from a Monte Carlo simulation using 
a parameterization ll38ll of the VGG model i39ll (see details in 
Ref. J24ll ). In each kinematic bin, the resulting systematic un- 
certainty is defined as the root-mean-square average of the five 
differences between the asymmetry amplitude extracted from 
the Monte Carlo data based on five GPD model variants l38ll 
and the corresponding model predictions calculated analytically 
at the mean kinematic values of that bin. 

Another source of systematic uncertainty comes from the 
relative shift of the squared missing mass distribution between 
e~ and e + data (see section |3). One quarter of the difference 
between the asymmetries extracted with standard and shifted 
missing-mass windows is assigned for this uncertainty. The 
background correction also makes a contribution to the uncer- 
tainty 12411 . There is an additional overall scale uncertainty aris- 
ing from the uncertainties in the measurement of the beam and 



target polarizations, which are given in Table Q] and stated in 
the captions of the figures and tables in the results section. Not 
included is any contribution due to additional QED vertices, as 
for the case of polarized target and polarized beam the most sig- 
nificant of these has been estimated to be negligible & The 
total systematic uncertainty in a kinematic bin is determined by 
adding quadratically all contributions to the systematic uncer- 
tainty for that bin. 

6. Results 

All of the asymmetry amplitudes in Eqs. [2"Tll2"7l are extracted 
simultaneously in a fit to the data. The results for beam-charge, 
charge-averaged or charge-difference single-spin asymmetry am- 
plitudes defined in Eqs. l2"Tll25l are compatible with those previ- 
ously published by HERMES l22l 12411 . They are not considered 
in this paper since these amplitudes are here extracted from a 
subset of previously analyzed data. 

The results for the Fourier amplitudes of the beam-charge- 
difference and charge-averaged double-spin asymmetries 
and y[™ +DVCS defined in Eqs.|26]and|27]are presented in Figs. [2] 
and[3]respectively, as a function of -f, xb, or Q 2 , while integrat- 
ing over the other variables. (In the HERMES acceptance, xb 
and Q 2 are highly correlated.) These values are also given in 
Tables |2 and [3] The 'overall' results in the left columns corre- 
spond to the entire HERMES kinematic acceptance. Figure [2] 
shows the leading amplitudes of the double-spin asymmetry re- 
lated to target transverse polarization combined with beam he- 
licity and beam charge, while Fig. [3] shows the amplitudes of 
the double-spin asymmetry, which relate to target transverse 
polarization and beam helicity only. The results for the various 
harmonics of the asymmetries and y[™ +DVCS were found 
to be compatible with zero within the total experimental uncer- 
tainties. The bottom row of each figure shows in each kinematic 
bin the estimated fractional contribution to the yield from asso- 
ciated BH production leading to a baryonic resonant final state. 
They are obtained from a Monte Carlo simulation using a gen- 
erator described in Ref. [24]. The two non-leading amplitudes 

A co S (0-fe)cos(2« ^ A sm(^0 s )sm(20) ^ ^ n = l m Eq ^ nQt 



LT,1 1 LT,I 

shown in Fig. [2] are found to be compatible with zero (see Ta- 
ble|2]l within the total experimental uncertainty. The correlation 
among all fitted asymmetry amplitudes is presented in Fig. |4] 

The curves in Figs. [2] and [3] represent results of theoretical 
calculations based on the GPD model described in Ref. 03911 . 
using the VGG computer program of Ref. |4in. A Regge ansatz 
for modeling the t dependence of GPDs 14211 is used in these cal- 
culations. The model J 39 |l is an implementation of the double- 
distribution concept yl |2J] where the kernel of the double dis- 
tribution contains a profile function that determines the depen- 
dence on controlled by a parameter b J43ll for each quark 
flavor. The theoretical calculations shown in these figures are 
obtained for the profile parameters b va i and fo sea equal to unity 
and infinity, respectively, which were shown to yield the best 
agreement with data for the beam-charge asymmetry ampli- 
tudes at HERMES Q. The leading amplitudes of the target- 



spin asymmetry extracted in Ref. Il24il have sensitivity to 
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the imaginary parts of the functions C l TP+ and C\ p _. The lat- 
ter has significant sensitivity to the CFF £, thereby providing a 
constraint on the total angular momentum of valence quarks ]44 , 
4511 . The width of the theoretical curves correspond to varia- 
tion of the total angular momentum /„ of w-quarks between 0.2 
and 0.6, with J c { = 0. In principle, the asymmetry amplitude 
A"j f ^ s ' sm ^ could provide a similar constraint through the real 
part of the function C\ p _, as can be seen from Eq. [10] Unfor- 
tunately, due to different kinematic prefactors, this amplitude is 
expected to be suppressed compared to those extracted from the 
asymmetry J^ m , and model calculations also indicate that it is 
much less sensitive to quark total angular momentum. 



7. Summary 

Double-spin asymmetries in exclusive electroproduction of 
real photons from a transversely polarized hydrogen target are 
measured for the first time with respect to target polarization 
combined with beam helicity and beam charge, and with re- 
spect to target polarization combined with beam helicity alone. 
The asymmetries arise from the interference between the deeply 
virtual Compton scattering and Bethe-Heitler processes. The 
asymmetries are observed in the exclusive region in missing 
mass that includes the proton together with baryonic resonances. 
The dependences of these asymmetries on —t, jcb, or Q 2 are 
investigated. The results for various harmonics of the asym- 
metries yi^ T and yi™ +DVCS were found to be compatible with 
zero within the total experimental uncertainties. Nevertheless, 
they may serve as additional constraints in global fits to ex- 
tract GPDs from measurements. The measured asymmetry am- 
plitudes are not incompatible with the predictions of the only 
available GPD-based calculation. 
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Figure 2: Charge-difference double-spin asymmetry amplitudes describing the dependence of the interference term on transverse target polarization in combination 
with beam helicity and beam charge extracted from hydrogen target data. These asymmetry amplitudes correspond to n = and n = 1 in Eq. [26] The error 
bars (bands at the bottom of the panels) represent the statistical (systematic) uncertainties. There is an additional overall 8.6% scale uncertainty arising from the 
uncertainties in the measurements of the beam and target polarizations. The curves show the results of theoretical calculations using the VGG double-distribution 
model f3^l4lll with a Regge ansatz for modeling the t dependence of GPDs 1421 . The widths of the curves represent the effect of varying the total angular momentum 
J u of u-quarks between 0.2 and 0.6, with Jj = 0. The bottom row shows the fractional contribution of associated BH production as obtained from a MC simulation. 




Figure 3: Charge-averaged double-spin asymmetry amplitudes describing the dependence of the sum of squared DVCS and BH terms on transverse target polar- 
ization in combination with beam helicity extracted from hydrogen target data. These asymmetry amplitudes correspond to n = and n = 1 in Eq. 1271 The error 
bars (bands at the bottom of the panels) represent the statistical (systematic) uncertainties. There is an additional overall 8.6% scale uncertainty arising from the 
uncertainties in the measurements of the beam and target polarizations. The curves and the bottom row of panels have the same meaning as in Fig. [2] 
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Table 2: Results for azimuthal Fourier amplitudes of the asymmetry 3v LT - An additional 8.6% scale uncertainty is present in the amplitudes due to the uncertainties 
of the beam and target polarization measurements. 
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Table 3: Results for azimuthal Fourier amplitudes of the asymmetry J?(™ + . An additional 8.6% scale uncertainty is present in the amplitudes due to the 
uncertainties of the beam and target polarization measurements. 



9 



A LT 


sin(<H s )sin(2<t>) 






A LT 


sin(c|>-ysin(<|)) 






A[_t 


cos^-c^cos^) 






A[_t 


cos((|)-4> s )cos(<|>) 


□ 


□ 


A LT 


cos(<M.) 






. BH+DVCS 

A LT sin(4>-t|) s )sin((|)) 






. BH+DVCS . . , » . 

A LT COS(0-0 s )COS((j)) 




□ 


.BH+DVCS , , , . 

A LT cos(<H a ) 






A 1 

M UT 


cos(4>-(|) s )sin(3i])) 






Aut 


cos(4>-<|> s )sin(2(M 






Aut 


cos((H s )sir#) 






Aut 


sin(4>-<|> s )cos(3(|>) 




□ 


A 1 


sin(4>-<|) s )cos(2i|)) 


□ 




Aut 


Bln^Jcoatt) 




□ 


Aut 


sin(<M>.) 


□ 




Alu 


sin(2<H 






Alu 


sin(<W 






.DVCS 
A UT 


cos(<j>-4> s )sin(2<|)) 


□ 




. DVCS 
A UT 


cos(<H s )sin(<M 




□ 


.DVCS 
A UT 


sin(i])-4i s )cos(2(()) 




□ 


.DVCS 
A UT 


sin^-^cos^) 


□ 




.DVCS 
A UT 


sin($-* s ) 




□ 


-DVCS 
A LU 


sin(<|)) 


□ 




A c 


cos(3<|>) 






A c 


cos(2<M 






A c 


cos(<|>) 




A c 


cos(0<M 







□ □ □ 





□ 




■ 




□ 




□ 












□ 












□ 


□ 










□ 






□ 




□ 


□ 


□ 


■ 


■ 




□ 


□ 


■ 


■ 




□ 


■ 






■ 


■ 


■ 


□ 


□ 


■ 


■ 


■ 


□ 


□ 


■ 


■ 


■ 


□ 


□ 



□ □ 
■ □ 



□ 
□ 



□ □ 

□ □ 
■ □ 



□ B ■ □ 

J I LlJ I L 



J L 





□ 






□ 






□ 


□ 


■ 


□ 


□ 


■ 




□ 


□ 


□ 


□ 


■ 




■ 


■ 


□ 


□ 


■ 




■ 


□ 


■ 


■ 


□ 




■ 


□ 


□ 


■ 


■ 




□ 


□ 


□ 


■ 






■ 


■ 


□ 


□ 


□ 






■ 


■ 


■ 


□ 


□ 










□ 







□ □ 

□ □ 



■ □ 

□ □ 



□ 



J I L 



J I I L 



□ 



□ 


□ 




□ 


□ 




□ 


■ 




□ 


□ 


□ 


■ 


□ 


■ 


□ 


□ 


□ 


■ 




■ 


□ 


□ 





J I I L 



<§- S- 



w tn ■=* 



S S. 

tn c *c c" 

8 JH jl £ ~k 

f 4- 4: -4 $ 



S S 5= ct" 



OOOO 



tn tn tn o 



£ £ £ £ o o 

tn tn tn tn o o 



tn o 
o o 

O tn 



§5 
< 



§5 
< 



< < < 



-5 -5 
< < 



-5 -5 
< < 



-6 -6 £5 
< < < 



< < < < < 



Figure 4: Correlation matrix for all fitted asymmetry amplitudes. The closed symbols represent positive values, while the open ones are for negative values. The 
area of the symbols represents the size of the correlation. 
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